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Abstract 

We calculate the rate for the deep-inelastic electroproduction of dijets at HERA. We 
study the weakening of the azimuthal (back-to-back) correlation between the jets, as x 
decreases, to see whether it can be used to identify BFKL dynamics from conventional 
fixed-order QCD effects. We show how this may give information on the transverse mo- 
mentum (for) dependence of the gluon distribution of the proton. 
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The production of dijets at the HERA electron-proton collider offers an excellent oppor- 
tunity to study the properties of the gluon distribution of the proton at small x. At lowest 
order, dijets^] are produced by the emission of a "hard" gluon from the initial or final state of 
the struck quark (the QCD Compton process jq —>■ gq) or by photon-gluon fusion (jg — > qq). 
The dijet events of particular interest are those in which the two jets tend to go in the virtual 
photon direction (in the jp centre-of-mass frame) but separated by a rapidity interval which is 
small relative to their large individual rapidities. In this configuration the proton-gluon fusion 
process dominates and small values of x g are sampled. Here x g is the longitudinal fraction of the 
proton's momentum carried by the interacting gluon. The process is shown diagrammatically 
in Fig. 1, where the dominant structure of the interacting gluon at small x g is exposed. We 
also show the transverse momenta Pit, Pit and kr of the outgoing (quark, antiquark) jets and 
the incoming gluon. 

Dijets can be produced at HERA via direct photons either by photoproduction (Q 2 ps 0) 
[|l], Q] or by deep-inelastic electroproduction (Q 2 of O(10 GeV 2 )). In the former case it is 
much more difficult to extract the direct photon events from the events in which the photon 
is resolved into its constituent partons 0. In fact it appears likely that a larger "clean" dijet 
event sample will be obtained for electroproduction, and so we study this process. We are 
especially interested in the properties of the gluon at small x g . However the values of x g that 
are sampled in any deep inelastic study exposing the final state are always greater than the 
Bjorken x. In our case x g ~ (1 + s/Q 2 )x, where V§ is the cm. energy of the produced dijet 
system. So the lower the values of pt, for which the jets can be clearly identified, the better. 

In principle, the calculation of the dijet production cross section requires integration over 
the complete phase space of the transverse momenta kr of the gluon. In the conventional 
approach we (i) restrict the integration over k\ to the region k\ <C p 2 T and (ii) let k\ = in 
the calculation of the subprocess cross section a^g — > qq), i.e. we evaluate a with the gluon 
on-mass-shell. As a result we get the familiar factorization formula, which symbolically is of 
the form 

x g9(x g , /i 2 ) <8> aiag -* qq, V 2 ) (1) 

where the scale fi 2 ~ p 2 T . In this way we can probe the conventional gluon distribution, g, of 
the proton. However, at small x g , the strong-ordering approximation is no longer applicable 
- we must keep the full k\ dependence of a and integrate over the full k\ phase space. As a 
consequence we must work in terms of the gluon distribution f(x g ,k^,fi 2 ) unintegrated over 
k\ @,@,||, that is 

2 ft* ak T 2 2 



Xg9(Xg,fJ. ) = / -T3-f(Xg,kT,H )• (2) 
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At small x g the function / satisfies the BFKL equation 0, 



df(x g , k\) _ 3a s 2 f dk'r 



d\og(l/x g ) 7i 2 J k' 2 



J{Xg,k T ) j{Xg,k T ) j{Xg,k 



+ 



\k%-k 2 T \ (4A# + A#)3 



(3) 



which effectively resums the large leading log(l/x ff ) contributions which arise from the sum 
of the gluon emission diagrams of the type shown in Fig. 1 together with the virtual correc- 
tions. Note that, at small x g , the gluon distribution / becomes independent of the scale /i 2 . 

: As is customary, we use the term "dijet" to refer to two jets produced in addition to the jet formed by the 
remnants of the proton. 
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Two characteristic features of the solution of the BFKL equation are (i) the leading small x g 
behaviour of the form 

f(x g , k%) ~ x- A (4) 

where A ~ 0.5, and (ii) a "diffusion" in k^ with decreasing x g which arises from the "random- 
walk" in the kx of the gluons as we proceed along the gluon chain || |J. The first property 
gives rise to a growth in dijet production with decreasing x g . However this behaviour can be 
mimicked by eq. (1) with a conventional gluon which evolves from a singular distribution at some 
starting scale, that is x g g ~ x~ x . The singular small x g behaviour is stable to the evolution. 
The second property is more unique to BFKL dynamics. The diffusion in kx manifests itself 
in a weakening of the back-to-back azimuthal correlation of the two outgoing jets of transverse 
momenta pyr and p 2 r- As x (and hence x g ) decrease, larger k^s are sampled and broader 
azimuthal distributions are expected. On the other hand the strong-ordering of the transverse 
momenta of conventional dijet production leads to a narrow distribution about the back-to- 
back jet configuration. Thus, in principle, a measurement of the azimuthal distribution offers 
a direct determination of the kr dependence of the gluon distribution f(x g ,k^). In practice 
the situation is not so clear. The azimuthal distribution will also be broadened by higher-order 
conventional QCD effects. To see whether these will mask the BFKL signal we therefore also 
compute the azimuthal distribution resulting from the emission of a third "hard" QCD jet. 

We begin by using BFKL dynamics at small x g to calculate the differential dijet cross section 
as a function of 0, the azimuthal angle between the transverse momenta pit and p2T of the 
two jets. That is we evaluate 



da Ana^ 



dxdQ 2 d<f) Q 4 x 



y\ dF T (x,Q\<P) , dF L (x,Q* 

, l ~y + tt) 7- + (i - y) 7— 
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(5) 



where, as usual, the deep inelastic variables Q 2 = —q 2 , x = Q 2 /2p.q and y = p.q/p.p e where p e , p 
and q are the four momenta of the incident electron, proton and virtual photon respectively, see 
Fig 1. The differential structure functions dFi/d(f> can be computed from the /c^-factorization 
formula, which is shown symbolically in Fig. 1. It is convenient to express the jet four-momenta 
in terms of Sudakov variables 

Pi = (1 - p)q' + aip + pit 

P2 = Pq +a 2 p + P2T (6) 

where q' = q + xp and p are basic lightlike momenta. Thus, since the jets are on-mass-shell, we 
have 

ai = { (ww ) *' a2 = {-Jo*-} x (7) 

where m q is the mass of the quark. The Bjorken variable x ~ Q 2 /2p.q'. The factorization 
formula for the differential structure functions is 

dri{X,Q ,(p) ^ f f dp lT dp 2T ., i2\n-Q(n 2 2 i\ /o\ 

Jl = 1^J d PJ — ^i f(x9,k T )F?[f3,p( T ,pi T ,(f)) (8) 

with i — T,L and where the "factorization" variables 

.2 ir& 



x g = x + a x + a 2 { > [1 + 4p iT /Q ]x j 

k 2 -, = p\ T + p\ T + 2p 1T P2TCOS(j). (9) 



The functions Tl which describe virtual photon-virtual gluon fusion, jg — >■ qq, at small x g are 
HOI 



n(P,A T ,pl T A) = e^ a ,f ( l-pf[^^-^ (11) 

where e q is the charge of the quark q and where the denominators 

A = p 2 T + m 2 + /3(l-/3)Q 2 . (12) 

The first two terms in {...} in (10) and (11) correspond to quark box contributions with pi 
being first a quark and then an antiquark jet, whereas the third term is the "crossed-box" 
interference term. Note that the apparent divergence of the integral in (8) at k? — (that is 
at (f> = it) is cancelled by the zeros of the functions / and Tt,li see (10) and (11). 

The only unknown in the determination of the differential cross section for deep-inelastic 
dijet production, da/dxdQ 2 d(j), is the gluon distribution f(x g , k?) that enters in (8). We calcu- 
late / by solving the differential form of the BFKL equation, (3), using knowledge of the gluon 
at x g = 10 -2 , as described in ref.0. The normalisation (though not the x g behaviour) of / is 
dependent on the treatment of the infrared k? region. However, the gluon distribution / can 
also be used to predict the behaviour of the structure function F 2 itself, via the inclusive form 
of the k<r- factorization formula (8), see ref. ||. Thus we can fix the infrared parameter in the 
determination of / so as to reproduce the low x measurements of F 2 at HERA |T2| . The result 



is shown by continuous curves in Fig.2, and corresponds to the choice k 2 a = 2 GeV 2 . We see 
that an excellent description of F 2 is obtained. Now that / has been fully specified in this way, 
we should be able to predict the azimuthal distribution reliably, provided, of course, that we 
sample sufficiently small values of x g for the BFKL solution to be appropriate. 

In Fig. 3(a) we show dFr/d(j) calculated from the /^-factorization formula, (8), for deep- 
inelastic dijet events with Q 2 = 10 GeV 2 for three different values of x. Each jet is required to 
have transverse energy squared, Ej>, greater than 10 GeV 2 . We notice a rapid increase in the 
dijet rate with decreasing x, and a weakening of the azimuthal back-to-back correlation. This 
last observation is more evident from Fig. 3(b) which shows the same distributions normalized 
to a common maximum value at <fi = n. The broadening with decreasing x is a manifestation of 
the diffusion in kx which is characteristic of BFKL dynamics. The detailed shape in the region 
of ~ it may not be reliable, since it corresponds to small values of the transverse momentum 
kr, of the gluon (see Fig. 1). Moreover hadronization effects will influence the distribution in 
this region. Rather we should study the normalised distribution, dFT/d<j), away from <ft = tc, 
say outside the interval 180 ± 20 degrees. 

The characteristic BFKL behaviour of the solution f(x, fcf-) of (3) only has a chance to set 
in for x < 10~ 3 . The precocious onset of this leading log(l/x) behaviour does indeed appear 
compatible with the striking rise of F 2 with decreasing x that has been observed at HERA, see 
Fig.2. However this is not conclusive evidence of BFKL dynamics. The F 2 data can equally well 
be described by Altarelli-Parisi (or GLAP) evolution. For example the dashed curves in Fig.2 
are the description obtained from a recent global structure function analysis based on (next-to- 
leading order) GLAP evolution from "singular" parton distributions at Q 2 = 4 GeV 2 JT^]. To 



distinguish BFKL dynamics from conventional QCD we must look into properties of the final 
state at small x, such as the weakening of the back-to-back correlations in the dijet events. 
However there is a price to pay. The BFKL dynamics is sampled at larger u x" for final state 
processes than is the case for the inclusive F2 measurement. In our example of dijets with 
£y(jet) > 10 GeV 2 and Q 2 = 10 GeV 2 , we see from (9) that we sample values of x g > 5x. 
Thus if f^Xg,^) assumes a characteristic BFKL behaviour for x g < 10~ 3 , then we anticipate 
that the broadening of the <fi distribution, with decreasing x, will only be relevant in the region 
x < few xlO -4 . This is near the limit of the region which is at present accessible at HERA. 
In Fig. 4 we therefore compare the dijet azimuthal distribution for x = 2 x 10~ 4 with that for 
x = 10~ 3 , that is two values of x which are appropriate for HERA. 

As mentioned above, the weakening of the back-to-back azimuthal correlations can also be 
obtained in a more conventional way from fixed-order QCD effects, in particular from 3+1 jet 
production. As usual the +1 refers to the jet associated with the remnants of the proton. Part 
of the 3+1 jet production is, of course, already included in the calculation based on BFKL 
dynamics, since the absence of strong-ordering in k? means that the gluonic ladder contains 
additional (gluon) jets. Before drawing final conclusions we must therefore compare our BFKL 
dijet predictions with the azimuthal distribution coming from conventional 3+1 jet production. 
To calculate the latter we use the PRO JET Monte Carlo ||14|| . We require two of the jets to 
have i?|(jet 1,2) > 10 GeV 2 and the third to have E^(]et 3) < 10 GeV 2 , chosen so that two 
jets are "visible" and the third (which may be either a gluon or a quark) is relatively "soft". 
The results are shown by the histograms superimposed on Fig.4. We have checked that the 
PRO JET predictions in the region \<f> — n\ > 20° are not sensitive to a reasonable variation of 
the cut-off, i/ij = Sij/W 2 > y , that is used to regulate the infrared singularities. 

We see from Fig.4 that the appearance of dijet events in the "tails" of the azimuthal distribu- 
tion, that is at angles such that \<f> — ir\ > 45°, at the predicted rate, will be a distinctive signal 
for BFKL dynamics. Nearer the back-to-back configuration the fixed-order QCD processes 
swamp the BFKL effect. Of course, since we work at the parton level and ignore the experi- 
mental problems of jet identification, only the BFKL signal can contribute for |0 — n\ > 60°. 
If the cut, _Ey(jet 3) < 10 GeV 2 , on the third jet is removed and the azimuthal distribution 
is plotted as a function of the angle, 0, between the two jets with the largest E^s then the 
PROJET prediction is essentially unchanged for \<f> — tt| ~ 20°, but is enhanced at larger angles 
with a steeper fall off towards the limiting angle |0 — n\ = 60°. 

The deep inelastic variables (x, Q 2 ) and jet cuts (E T > E 0T ) have been chosen in an attempt 
to optimize both the event rate and the BFKL signal at HERA. Clearly if we were able to go 
to smaller x (or, to be precise, smaller x g ) then the BFKL effect would be more pronounced, 
see Fig. 3. Now the observable deep inelastic region at HERA lies in the domain Q 2 /x < 10 5 
GeV 2 . On the other hand, we see from (9) that x g > x + 4(x/Q 2 )Eq T , when two jets with 
E T > E 0T are recorded. A higher jet threshold, E 0T , has the advantage that the fixed-order 
QCD contribution is suppressed relative to the BFKL signal, but (since Q 2 /x is bounded) then 
a higher x g is sampled. The lower plot in Fig.4 shows the results for Q 2 = Eq T = 25 GeV 2 and 
x = 5 x 10~ 4 , for which x g is increased by a factor of 2.5 in comparison to that sampled for 
Q 2 = E 2 T = 10 GeV 2 and x = 2 x 10" 4 . 

To sum up, the main aim of our paper is to quantify the observation that BFKL dynamics 
weakens the azimuthal correlation of the q, q jets produced in small x deep-inelastic scattering 
via the photon-gluon fusion mechanism. We are able to obtain an absolute prediction since 



the parameter which specifies the infrared contribution to the BFKL equation is chosen such 
that the measurements of F 2 at HERA are reproduced. (That a physically reasonable choice 
of a single infrared parameter suffices to describe the observed small x behaviour of F 2 is a 
far from trivial test of BFKL dynamics.) For the dijet events we find a substantial broadening 
of the azimuthal distribution and an increase of this broadening with decreasing x. However, 
at HERA energies, we find that the fixed-order QCD contribution from 3+1 jet production 
exceeds the BFKL signal near the back-to-back configuration. Nevertheless at sufficiently large 
values of |0 — 7r| BFKL dynamics dominates and give rise to a distinctive "tail" to the azimuthal 
distribution at an observable rate. In this way dijet production at HERA offers an opportunity 
to study the kr dependence of the gluon distribution of the proton. 
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Figure Captions 

Fig. 1 A diagrammatic representation of dijet production by photon-gluon fusion, jg — > qq, at 
small x g . The function / is the (unintegrated) gluon distribution of the proton. The 
cross section is given by the kr- factorization formula (8), which has the symbolic form 
o = f <8> F ■ T denotes the quark box (and crossed box) contribution. 



Fig. 2 The measurements of Fi at HERA (preliminary data from the 1993 run]T2[) shown to- 
gether with the BFKL description || (continuous curves) and the MRS (A) parton analysis 
fit |jT~3|1 (dashed curves). The measurements of the HI collaboration at Q 2 = 65 GeV 2 are 
shown on the Q 2 = 60 GeV 2 plot. 

Fig. 3 (a) The distribution dFx/dcj) predicted by the BFKL /^-factorization formula, (8), for 
deep-inelastic dijet events with Q 2 = 10 GeV 2 and £f.(jet) > 10 GeV 2 . (b) The distribu- 
tions normalized to a common maximum. 

Fig. 4 The curves are the differential cross section for dijet production predicted by BFKL dy- 
namics, (5)-(12), whereas the histograms correspond to 3+1 jet production as determined 
by the PRO JET Monte Carlo using MRS(A) partons fl3| . In the first case, the broaden- 
ing of the azimuthal distribution arises from the fcr dependence of the gluon distribution 
f(xg,kj<) found by solving the BFKL equation, (3). For 3+1 jet production we assume 
that the third jet has E\ < 10(25) GeV 2 in the upper (lower) plot. 
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